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Q^i Abstract 

C^ I The phase space of general relativity is extended to a Poisson manifold by in- 

\^ • elusion of the determinant of the metric and conjugate momentum as additional 

independent variables. As a result, the action and the constraints take a poly- 
nomial form. We propose a new expression for the generating functional for the 
\j^ ' Green's functions. We show that the Dirac bracket defines a degenerate Poisson 

^^ ■ structure on a manifold and the second class constraints are the Casimir functions 

.^ ! with respect to this structure. As an application of the new variables, we consider 

f^ ' the Priedmann universe. 

o 
^; 1 Introduction 

{3JQ, Canonically formulating any model of mathematical physics is the most important step 

when analyzing equations of motion, in particular, when setting and analyzing the Cauchy 
problem. It also provides the basis for canonically quantizing models. Canonically for- 
mulating general relativity is technically involved, and many papers are devoted to this 
problem. We mention only a few of them. Dirac first formulated general relativity self- 
consistently in the second-order formulism pT. He took the metric components Qap as 
the independent variables and showed that the Hamiltonian of the gravitational field is 
equal to a linear combination of constraints. Afterwards, Arnowitt, Deser, and Misner in 
the series of papers resulting in the review [2j essentially simplified the calculations and 
clarified the geometrical meaning of the canonical momenta, expressing them in terms of 
the extrinsic curvature of a spacelike hypersurface imbedded in a four-dimensional space- 
time. The expression for the Hamiltonian was found in the first order formalism when the 
metric g^p and the symmetric affine connection T^a/3}^ considered independent variables. 
In essence, this approach simplified the calculations. In the pioneering papers pj |2] , the 
constraint algebra was not calculated explicitly, but the constraints were shown to be 
consistent with the equations of motion (the first class constraints). In j2], the role of 
boundary terms was also analyzed, and the total energy of the gravitational field for an 
asymptotically flat space-time was defined in terms of the surface integral. The role of 
boundary terms was discussed in a more general form in |3|. 
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In j3] , general relativity was canonically formulated and a generalization of the Schrodinger 
equation for the wave function of the universe, which later was called the Wheeler-DeWitt 
equation, was considered in detail. The constraint algebra in general relativity was first 
calculated explicitly there. 

The constraint algebra for any model invariant with respect to general coordinate 
transformations was obtained in ^ assuming that the model is self-contained and that 
the constraints generate the general coordinate transformations for the canonical variables 
(see also [?]). We note that the self-consistency of a model (the closedness of the con- 
straint algebra) is a very strong assumption. Conversely, for a given model, the constraint 
algebra, because it is not known beforehand, must be calculated explicitly to prove the 
self-consistency of the model. 

Dirac [7] and Schwinger jS] started the investigation of the vielbein Hamiltonian for- 
mulation of general relativity using the time gauge for simplicity. The Hamiltonian for- 
mulation in a general case without gauge fixing was given much later because of serious 
technical difficulties piTU]. 

The Hamiltonian formulation of general relativity contains constraints that are non- 
polynomial on space-section metric components, and this is an essential obstacle to an- 
alyzing and quantizing the theory of gravity. The polynomial Hamiltonian formulation 
given by Ashtekar ^j attracted much interest in recent years. He proposed using complex 
variables in the extended phase space, which are tensor densities and lead to polynomial 
constraints. Here, we consider a different extension of the phase space ^2] with the metric 
determinant and its conjugate momentum considered additional variables. We show that 
the Poisson structure on the extended phase space is degenerate and that the initial phase 
space is mapped on a subspace of the extended phase space by the canonical transfor- 
mation. All new canonical variables are real tensor densities, and the constraints take a 
polynomial form. 

We propose a functional integral over a Poisson manifold as a new expression for the 
generating functional for the Green's functions. This form of the integral reduces to the 
standard expression for the generating functional over the phase space |13^ after integra- 
tion over the additional variables, which is removed by two supplementary (5-functions. 
We prove that the corresponding Jacobian of coordinate transformation is equal to unity. 

The plan of the paper is as follows. In Sees. |2Hni we describe the transition from the 
Hilbert-Einstein Lagrangian to the Hamiltonian in detail to avoid sending the reader to 
the original papers, where a greater part of the calculations is usually omitted. Moreover, 
we consider a general case of affine space-time geometry in Sec. El when describing the 
geometry of hypersurfaces. In this case, the antisymmetric part of external curvature 
of a hypersurface is defined by the torsion tensor. This is important for the canonical 
formulation of general relativity in the vielbein formulation and for the Hamiltonian for- 
mulation of models with absolute parallelism. The canonical transformation between the 
phase space of general relativity and the submanifold of the extended Poisson manifold 
is described in Sec. |7| We show that all constraints and the action of the model take a 
polynomial form in the extended space, and we compute the constraint algebra. We pro- 
pose the expression for the generating functional for the Green's functions on the Poisson 
manifold in Sec. |H1 In Sec. IHlas an application of the new variables, we consider the case 
of a homogeneous and isotropic universe, a case where all the calculations can be easily 
checked. 



2 ADM parameterization of a metric 

To analyze the Hamiltonian structure of the general relativity equations Arnowitt, Deser, 
and Misner used the special parameterization of the metric (ADM-parameterization), 
which essentially simplifies calculations 2J. We consider a manifold M, dimM = n 
equipped with a metric of Lorentzian signature (H — ...—). We deliberately do not re- 
strict ourselves to the most important case of four- dimensional space-time because gravity 
models in higher and lower number of dimensions have attracted much interest recently. 
Let {x"}, q; = 0, l,...,n — 1 denote the local coordinates. We choose the time coordinate 
t = x", and then {x"} = {x", x'^}, /i = l,...,n — 1. In what follows, the letters from the 
beginning of the Greek alphabet (a, /5, . . . ) range all index values, while the letters from 
the middle (/x, z/, . . . ) range only the space-related values. This rule is easily remembered 
from the inclusions {/i, v, ■ ■ ■} C {a, /3, . . . } and {1,2,...} C {0, 1,2,...}. The ADM 
parameterization of the metric has the form 

where g^i, is the metric on (n— l)-dimensional manifold sections x° = const. In the chosen 
parameterization, we introduced the same number of functions N and A^^ instead of the 
n metric components containing at least one time index qqq and qq^. Here, A^'' = g''^N^, 
where g''^ is the (n — 1) x (n — l)-matrix inverse to (7^,^: 

which we call the inverse metric on the sections x° = const. In what follows, we always 
raise the space indices using the inverse metric g'''^, which is marked with the hat and does 
not coincide with the space part of the metric g"'^ inverse to gap, g''^ 7^ g''^- The function 
A^ = A^(x) is called the lapse function, and the functions A^ = A^(x) are shift functions. 
Without loss of generality, we assume that the lapse function is positive (A^ > 0). In this 
case, the ADM parameterization of the metric (0) is in one-to-one correspondence. The 
interval corresponding to the parameterization (P) has the form 

ds^ = N^df + g^^idx'' + N''dt){dx'' + N^dt). 

We assume that the coordinate x*^ = t is the time, i.e., that the vector c^o tangent to 
the coordinate x" is timelike. Formally, this condition is written as 

{do,do) = goo = N^ + N''N^>0. (2) 

In this case, the metric gai3 has the Lorentzian signature if and only if the matrix 

A-^A"^ 
^^^ ~ m + NPN, ^^^ 

is negative definite. We note that the metric g^u induced on sections x*^ = const may not 
be negative definite. This means that sections x° = const are not spacelike in a general 
case. In what follows, we additionally assume that the coordinates are chosen such that 
all sections x*^ = const are spacelike, i.e., the metric gf^i, is also negative definite. This is 
convenient for posing the Cauchy problem when the initial data are given on a spacelike 
surface and we consider their evolution in time. 



Similarly, a metric on a Riemannian manifold can be parameterized with a positive- 
definite metric. For this, it suffices to explicitly choose any coordinate instead of time. 
The metric inverse to (0) is 



nl3 







(4) 



The space matrix in the lower right block 
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flU 



9^" + 






(5) 



is inverse to metric (jH)), as can be easily verified. This means that the negative definiteness 
of metric Q is equivalent to the negative definiteness of the matrix g'^'^. 

We note that if the metric on a manifold M 
has the Lorentzian signature, then the condi- 
tion that all sections x° = const are spacelike 
is equivalent to the condition A^^ > 0. Indeed, 
the negative definiteness of the inverse matrix 
gtj-i' follows from that of Qf^^. Then the negative 
definiteness of the matrix 
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fii/ 






follows from Eq. (0). In turn, this is equivalent 
to the condition 5f°° > or A^^ > 0. 

We consider a simple example to show de- 
tails that may arise for ADM parameterization of a metric. 

Example. We consider the two-dimensional 

Minkowskian space-time M^'"^ with the Cartesian coordinates t, x. We introduce the new 
coordinate system ^, rj depending on two real parameters a and h (see Figure) 

^ = t + ax, ri = t — bx, |a| 7^ 1, |6| 7^ 1, a + b^O. 

We can easily obtain the formulas for the inverse transformation 

h^ + arj i — V 



t 



a + b 



X 



a + b 



The metric has the form 

ds'^ = dt^ — dx^ - 



(a + by 



[(62 - l)rf^2 ^ 2{ab + l)did7] + (a^ - i)drf] . 



in the new coordinates. 

We now analyze the ADM parameterization of the metric in the coordinates x^ 
$,, x^ = r]-. 



9oo 



(a + 6)2' 
The lapse and shift functions are 



N' 



9oi 



ab + 1 



Ni 



9ii 



{a + by 



ab + 1 
(a + 6)2- 



The inequalities |6| > 1 and |a| < 1 follow from the respective conditions goo > and 
Qii < 0. We see that these conditions are necessary and sufficient for the coordinate lines 
^ and Tj to be respectively timelike and spacelike. It is easy to verify the equivalence of 
the conditions 

NiNi 1 

Using the formula for the determinant of block matrices, we obtain the expression for 
the determinant of metric (^Q) : 

detgap = N'^detg^,^. (6) 

We hence have the expression for the volume element 



e = Ne, e= ^J\detgaf^\, e = y | det gf^^ | . (7) 

This formula is a generalization of the well-known school rule: the volume of a prism is 
equal to the product of the base area and the height. In this case, the base area is e and 
the height is the lapse function A^. 

The following formulas, which can be varified straightforwardly, are useful for calcu- 
lations: 



g'^'Qu. = <5^ + 






g^ Qi^u =n-l + 



3 Geometry of hypersurfaces 

In the Hamiltonian formulation of gravity models, we take the space-time as a family 
of spacelike hypersurfaces x^ = const parameterized by time. In other words, in each 
instant, the space is a hypersurface embedded in the space-time. It is useful to know what 
geometry arises on spacelike hypersurfaces because equations of gravity models define 
geometry of the whole space-time. In the present section, we approach this problem 
from a general standpoint, assuming that an arbitrary affine geometry is given on the 
embedding manifold without assuming that the metric signature is Lorentzian. 

We consider an (n — l)-dimensional hypersurface U embedded in a n-dimensional 
manifold M: 

/ : U ^ M. (9) 

We let x", a = 0, 1, . . . , n — 1, and u\ i = 1, . . . ,n — 1, denote the respective coordinates 
on M and U. Then the embedding of U in M is locally given by n functions x"(-u), 
which we assume to be sufficiently smooth. An arbitrary vector field {X*} G T(U) on the 
hypersurface is mapped on the vector field {X"} G T(M) on M by the map differential 

/, : X = X'di G T(U) ^ Y = r"9„ G T(M), 



where 

The Jacobi matrix e"j of the transformation / is rectangular of the size n x (n — 1), has 
the rank n — 1, and is obviously irreversible. It is defined not on the whole manifold M 
but only on the hypersurface U. In addition, we note that the Jacobi matrix is a vector 
and covector with respect to the respective coordinates transformations on M and U. The 
pullback of the map / maps each covector field on the image /(M) C M in the covector 
field on U: 

/* : A = rfxM„GT*(M) -^ B = du'BieTi\]), 

where 

R- — 4 p" 

We identify the hypersurface U with its image U = /(U) C M in what follows. 

The 1-form n = dx'^na defined on the hypersurface U by the system of algebraic 
equations 

nae°'i = 0, i = l,...,n — 1, (10) 

specifies the field of (n — l)-dimensional subspaces tangent to U in the tangent bundle 
T(M). These equations have a unique solution up to multiplication on an arbitrary 
nonzero function because the rank of the Jacobi matrix is equal to n — 1 as a consequence 
of the definition of embedding. 

The Jacobi matrix e"j defines a set of ri — 1 vectors Cj = e'^ida in the tangent spaces 
T2.(M), X G U; these vectors form the basis of the space tangent to the hypersurface. 

This is all we can say about a hypersurface U if there is only embedding Q. The 
theory becomes much richer in content if there are additional structures on M. We discuss 
this question in detail. 

Let the afiine geometry be given on M, i.e., a metric gap and an affine connection 
Ta/^^- We consider what geometry arises on the hypersurface U C M. The pullback of 
the map /* induces a unique metric on the hypersurface: 

/* : 9af3 -^ 9ij = 9ai3e"'ie'^j. (11) 

The existence of the respective metrics Qais and gij on M and U allows lowering and raising 
the indices of the Jacobi matrix: 

ea = 9af3e^j9'^, 

where g^^ is the metric inverse to gij. This matrix projects an arbitrary vector from T^.(M), 
X G U, into the space tangent to the hypersurface T(U) 

We now define the connection on the hypersurface U C M by the relation 

V,X' = (V,X^)e^e/. (12) 

Opening this relation leads to the expression for the induced connection on the hypersur- 
face U in the coordinate form: 

f,/ = id^x^' + r,^^e^e^,)e/. (13) 



This connection is unique. We note that if the original connection F^./j''' is symmetric, 
then the induced connection r^^ is also symmetric. As a consequence of Eq. ()13|) . the 
torsion tensor T^/j''' = Tap'^ — Tj3a^ on M induces the torsion on the hypersurface 

Furthermore, the connection on U is uniquely defined only in the case where the metric is 
given on M in addition to the connection. All geometric objects related to the hypersurface 
and constructed using only the induced metric gij and connection Fjj'^ are marked by the 
hat in what follows. 

The metric Qajs and connection Fq,^''' on M define the unique metric Qij and connection 
Tij'' on the hypersurface U C M. The converse statement is not true. If a metric and a 
connection are given on the hypersurface U, then they do not induce the geometry on M 
uniquely. This is clear because the dimension of the hypersurface is less than that of the 
many fold itself. 

Straightforward calculations yield the expression for the covariant derivative of the 
induced metric on the hypersurface: 

^igjk = diQij - Ti/gik - TikQji = {Vagi3-y)e°'ie'^je^k- 
This relation gives the expression for the nonmetricity tensor on the hypersurface 

^ijk ^aP^(^ jC j'C fc. 

In particular, if the connection Tai3^ on M is metrical [Qap-y = 0), then the induced 
connection Fjj'^ on U is also metrical {Qijk = 0). 

The existence of metric Qais allows forming the unit vector field n = n°'da orthogonal 
to the hypersurface. It was already noted that the system of equations riae'^i = defines 
the 1-form dx°'na up to multiplication on an arbitrary scalar function. We use this 
arbitrariness for the vector ra" = g°'^np to have the unit length in every point: (ra, n) = 
n°'n^ga/3 = 1- This vector is orthogonal to all vectors tangent to the hypersurface by 
construction: 

{n, d) = n"'e'^igaf3 = n^e'^i = 0. (15) 

If the hypersurface is given on a manifold, then there is a natural basis {n, Cj} in the 
tangent space T(M) defined by this hypersurface. This basis is defined only in points of 
the hypersurface but not on the whole manifold. The dual basis {n = dx°'na, e* = dx°'ea^} 
in the cotangent space T*(M) corresponds to it. Then an arbitrary vector X and a 1-form 
A can be decomposed with respect to this basis: 

X° = X^n" + X'e'^i, X^ = X°n«, X' = X"e«\ 

Aa = A±na + AiCa , A± = AaTl , Ai = A^e i. 

A tensor of an arbitrary rank can be decomposed similarly. In particular, a covariant 
second rank tensor has the decomposition 

Aai3 = A_ij_nanp + Aj_inae'i3 + Ai^e\np + Ajj-e^e^/j, 

where 

^±± = A^pn^'n'^, Au = A^pn^'e^i, Ai^ = Aape^iU^, Aij = A^pe'^iC^p 



We can easily verify that the decomposition of the metric is essentially simpler: 

9ap = n^np + CaCp^gij. (16) 

A similar decomposition holds for the inverse metric: 

The summation over Latin indices for the Jacobi matrix follows from the definition of the 
inverse metric g'^^gp-y = S^: 

e"^e/ = 5| - n"n^. (18) 

As a consequence of Eq. ()15|) and the definition of the inverse induced metric g^-'gjk = 
61, we have the equality 

e'^iej = 51 (19) 

where the summation is over the Greek indices. Using this rule we obtain the represen- 
tation for the inverse induced metric 

which follows from ()17|) . Metric p6|) and its inverse ()17p in the basis n, Cj have the block 
diagonal form: 

'1 0\ A 

,0 g^,)' \Q g^\ 

This allows raising and lowering the corresponding indices. For example, if Xq = X^^gpa, 
then Xx = X^ and Xi = X^gji. 

Induced metric (fTT|) and connection ()13p define the internal geometry of the hypersur- 
face U C M. In particular, the induced connection yields the internal curvature tensor of 
the hypersurface: 

Rijk (r) = diTjk — Tik'^Tjm —(«<-> j). 

The embedding / of the hypersurface allows defining an additional important object 
which is called the external curvature of the hypersurface, 

K,, = -V^^n^e^e^,-, (20) 

which is equal to the covariant derivative of the normal projected on the space tangent to 
the hypersurface up to a sign. In contrast to the internal curvature tensor, the external 
curvature is a second-rank tensor, which has no symmetry in the indices in the general 
case. This tensor characterizes the variation of the normal when it is translated parallel 
along a curve on the hypersurface. Expanding this definition and using (jlU|) we obtain 

i^^i = n„(9jx" + r^/eV,)- 

The antisymmetric part of the external curvature tensor is given by the torsion tensor, 

K,, - Kji = 2Kiij^ = n<,T^/e^e^, = T,/. (21) 

As a consequence, the external curvature is symmetrical if and only if the connection T^^" 
has no torsion. 

The covariant derivative of the Jacobi matrix is 

V.e",- = eMdfse'^, + T^/e^,) - f,/e"fc, (22) 



where we use the connection F^/?'^ on the whole manifold along with the connection f jj*^ 
on the hypersurface. Simple calculations show that this covariant derivative has only the 
normal component and is proportional to the external curvature: 

Vie%- = n^'Kij. (23) 

This relation is known as the Gauss-Weingarten formula. As a consequence, we have one 
more representation for the external curvature tensor: 

Kij = n^Vie'^j. (24) 

The full curvature tensor Raises of a manifold M projected on the hypersurface can be 
expressed in terms of the internal curvature tensor Rijki constructed for induced metric 
fill!) and connection (fTSj) and the external curvature tensor. For this, we consider the 
commutator of covariant derivatives of a vector field, which is given by the curvature and 
torsion tensors, 

[Va, V/3]X.y = —Ra/S-y Xs — Ta/3 V^X^ = 

= —Rap-ylX — Raf3y±X — Tafj V^X^ — Tafj V^X^, 



where we first compute the covariant derivatives in the right-hand side and then project 
them on the hypersurface: V^X^ = e"/VaX^, V_lX^ = rfVaX^. To project this relation 
on the hypersurface, we first project the covariant derivative: 

WiX, = e^(V,X;3)e^- = V.X,- - X^Ki^, 

where 

VjXj = diXj — Tij Xk 

is (n — l)-dimensional covariant derivative on the hypersurface. The second covariant 
derivative is projected similarly: 

VjVjXfc = VjVjXfc — \/±XkKij — \/jX±Kik = 

= ^NjXk - ViX^Kjk - VjX^Kik - XjViKjk - 'V±XkKij, 

where 

V,ir^ = e^(V«X;5)n^ = V,X^ + X^K,,. 

The antisymmetrization of the obtained expression in the indices i,j yields the projection 
of commutator (|25p on the hypersurface: 

[Vj, Vj]Xfc = —RijkiX — Rijk±X — Tij V/Xfc — Tij V^X^. 

Taking the independence of X' and X^ and expressions (fT^ and (PT|) for the torsion and 
curvature tensors into account, we obtain the expressions for the projections of the full 
curvature tensor on the hypersurface: 

Rijki = Rijki + KikKji — KjkKii, (26) 

Rijk± = "^iKjk — S/ jKik + Tij Kik- (27) 

The obtained relations are generalizations of the Gauss-Peterson-Codazzi equations to 
the case where an arbitrary affine geometry with nonzero torsion and nonmetricity is 
given on the embedding manifold M instead of a Riemannian geometry. 
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To conclude this section, we compute the normal components G±± and G±i of the 
Einstein tensor 

GafS = Rap — -^QapR- 

in the Riemannian geometry where the torsion and nonmetricity equal zero. First, we 
compute the scalar curvature 

R = g'^^g^'R^p^s = 2R±^ + g''g''R^Jkl, 

where -R_l^ = g^-'Ri±j± is the normal component of the Ricci tensor and the expression for 
the inverse metric (fTTjl is used. As a consequence of the Gauss-Peterson-Codazzi ()26|) . 
we obtain 

g"'g''R,,M = R + K^-K''K,,, 

where R is the scalar internal curvature of the hypersurface and K = g^^Kij is the scalar 
external curvature of the hypersurface. As a consequence, we obtain the expressions for 
the normal components of the Einstein tensor 

2^ '^' (28) 

G±i = WjK^ - WiK. 

It is important that these components of the Einstein tensor do not contain derivatives 
normal to the hypersurface V_l of the induced metric and external curvature tensor at all. 
In the Hamiltonian language, this means that time derivatives are absent and that Ein- 
stein's equations G^±_ = and G±i = represent the constraints because the components 
of the external curvature K^^ are shown to be proportional to the momenta canonically 
conjugate to the induced metric gij in the next section. 

4 Curvature in the ADM-parameterization of the met- 
ric 

The ADM parameterization of metric (0) is convenient for the canonical formulation of 
general relativity in which the metric components gap and canonically conjugate momenta 
p°'^ are the independent variables. The passage from the Lagrangian to the Hamiltonian 
needs a relatively tedious calculations, which are given here. 

To essentially simplify calculations, we use the results of the preceding section. Namely, 
the sections x° = const of the space-time M yield the family of the hypersurfaces U C M 
which are spacelike by the assumption. We take the space coordinates as the coordinates 
on the hypersurfaces 

As a consequence, we loose the freedom of independent coordinate transformations on the 
space-time M and on the spacelike hypersurface U, but many formulas became simpler. 
The Jacobi matrix of the hypersurface embedding in the case under consideration is 

where 0,^ denotes the row consisting of n — 1 zeroes. The embedding induces the metric 
gfj^i, on the hypersurfaces according to formula (jllj) . 
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We now construct the vector field n = n^da orthogonal to the family of hypersurfaces 
x^ = const. As a consequence of orthogonality condition {n,X) = 0, where X = X^d^ is 
a vector tangent to the section x^ = 0, we obtain 

Moreover, if we set n° = 1/A^, then the length of the normal vector is equal to unity 
{n? = 1). The unit vector orthogonal to a section a;° = const has the form 

n = ^{do-N^d,) (29) 

and is always timelike. The corresponding orthonormal 1-form is 

n = dx^N. (30) 

An arbitrary vector on M can be decomposed with respect to the basis {n, e^}. In 
particular, we have the decompositions for vectors and 1-forms 

X" = X^n" + X'^e^", X^ = X^n^ + X^e^„, 

where 

x^ = x°x, x^ = x°x^ + x^. 



1 



Xi. = -(Xo - X^X^), x^ = x^. 



The representations for the metric (|16|) of the whole space-time and its inverse (|17|) 
are 

gap ^ ^a^/3 ^ ^M^g^ag^/3_ ( ) 

Connection (fT5|) induced on the hypersurfaces is the Christoffel symbols T^^p computed 
for the space metric g^u. 

In the ADM-parameterization of the metric, external curvature tensor (j^UI) on a hy- 
persurface x^ = const has the form 

K^u = r^.°x = ^(v^iv, + v,x^ - V), (32) 

where the dot denotes the differentiation with respect to time, g^u = dog^i, and V^X^ = 
dfj.Ni, — Vfju^Np. The external curvature tensor is symmetrical (X^j, = Xj,^) because the 
torsion vanishes in the metric formulation of general relativity. In what follows, we need 
the trace of the external curvature tensor 

X = X/ = ^'^^X^,. 

All time derivatives of the space part of the metric g^ji, are conveniently expressed in 
terms of X^j, when computing the curvature tensor Ra/S'yS of the space-time M. Moreover, 
to exclude the second time derivatives (jfj,,^, we need the time derivative of the external 
curvature tensor 



K = 

^'^ 2X 
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X 
K 



where 



v^iv, = d^N, - r^/iv. 



and where we must substitute the expression for g^^^ in terms of K^i,. 

We now compute the curvature tensor Rap-yS for the metric of form ((H) . Straightforward 
calculations yield the linearly independent Christoffel symbols 



Too" = ^ (^ + ^'^p^ + N^N'^K^ 

Too'' = g^^iN^ - Nd^N - N^V^Np) - — (n + N^dpN + N^N'^K^ 






pa j 1 






r ° = —K 

r ^ = r p K 

In what follows, we need the following combinations of the Christoffel symbols 



r« = r,/, ^^^^r^/ 



Simple calculations give 



N 
N 



To = - + V^iV^ - NK, 



r =f +M 

1 _ 1 



+ -:jg'"'{Np - NdpN - N'^VpN, - 2N'^V.Np + 2NN'^Kp^) 



1 
We also write the formulas for the time derivatives of the Christoffel symbols 



(33) 



g^^r^,' = -K + —{N-mdpN), (34) 

/^r,/ = l^r + -j^) Tp.- -—K- —(AT - NPdpN) + 



do^fiup — -iy pCjup + ^uQpp — ^pQpu) + ^pu'^Qpa, 

5of ^." = \r''iy^.gup + V.gpp - Vpgpu). (35) 

dotp = Ig'-'Vpg^p. 

The time derivatives g^^u are also eliminated from these expressions using relation (jH^ . 
We now compute the linearly independent components of the curvature tensor: 

Ropou = - Nkp, + R^p^^NPN'' + NNP{VpK,p + V.Kpp - VpK^,) + 
+ NV^W.N + Kp^N^N'^Kp^ + N{KpPt,Np + K/VpNp)- 
- N'KpPKp - NPN^KppK,, (36) 

Rpupo =R^.up.N'' + N{VpK,p - V.Kpp) + {KppK,, - K,pKp^)N'' 

R =R + K K - K K 
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where we use the formula for the commutator of covariant derivatives 

The components of the curvature tensor having at least one time index seem simpler 
when expressed in the basis n, e^: 

(37) 



-'^flUpl. ^ fjb^Up ^ V^ Ijbp- 



The components of the curvature tensor R^ypa and Rpup± were actually obtained in the 
preceeding section ^I^ . (j77|) without straightforward calculations. 
The Ricci tensor has the linearly independent components 

Roo = - N^^Kp, + Rp.N^'N- + NN^{2V,K\ - d^K) + NV^VpN+ 









Rop ~ [-Kp, + N'^V.Kp, + V^V.iV + ir/V.iV, + ir/V^iV, 



Rp. =Rpu + ^ -i^M- + V^V.iV + i^/ V.iV, + K/V^Np ] + 



+ f^^.iV' + iV (V.K% - 9^k) + +Kp,N''K - 2K/K,piV^ 

1 

iv 

+ ^VpKp, + Kp^K - 2KpPK,p. 
For reference, we also write the Ricci tensor components with respect to the basis n, e^: 

R±p =^uK\ - dpK. 
Finally, we compute the scalar curvature: 

R = R+^ [-r^'Kp, + V^VpN + 2K^'''VpN, + N^dpK^ - SK'^'^Kp, + K^. (40) 

5 The Hamiltonian 

The scalar curvature contains second derivatives of the metric components with respect 
to time as well as the space coordinates and is therefore not suitable for canonically 
formulating general relativity. It suffices to eliminate only the second time derivatives 
from the Lagrangian to obtain a canonical formulation. The Lagrangian density takes the 
simplest form after adding the boundary term: 

C^^u = NeR + 29o (eK) - 2dp (eg'^'d^N) . (41) 

Straightforward calculations yield the expression 

Cj^uM = Ne (^K^'Kp, -K^ + ny (42) 
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The transition to the Hamiltonian formahsm is now easy. First, the ADM Lagrangian 
does not contain time derivatives of the lapse and shift functions A^ and A*"^. This means 
that the theory contains n primary constraints 

their number coinciding with the number of independent functions that parametrize dif- 
feomorphisms. 

The momenta canonically conjugate to the space metric g^i, are proportional to the 
external curvature tensor, 

We note that the momenta are not tensors with respect to the coordinate transformations 
x'^ but tensor densities of the degree —1, as is the determinant of the vielbein, which degree 
is 

dege = —1 

by definition. 

To eliminate the velocities g^iy from the ADM Lagrangian, we decompose the momenta 
on the irreducible components, extracting the trace from p'^'^ 

p/^^=pM^ + ^^j9^M^, (45) 

n — 1 

where we introduce the trace of the momenta 

P = P^'g^.. = e(n - 2)K (46) 

and the symmetric traceless part 

We can now solve Eq. (jl^ for the velocities using relation ()32p . 

2N f 1 \ -^ 

V = — [Pt^u- / _ 2) P3f^f 1 + "^ti^u + ^uN^. 

Simple calculations yield the Hamiltonian density 

^ = P^^'g^.u - -Cadm = NH^ + N^H^ + 2d^{p^''N,), (47) 

where 

H^ = \ (p^^'p,. - Y^y\ - eR, 

e V {n-2) J (48) 

H^ = -2V,p\ = -2d,p\ + d^g.pp'^P, 

and Pfj,u = gfj.pgiyaP'"^ ■ We note that the covariant derivative of the momenta contains 
only one term with Christoffel symbols because the momenta are tensor densities. The 
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expressions for H^ and H^ are proportional to the G^±_ and G^^ components of Einstein 
tensor ()28|) . which justifies the chosen notations. 

Dropping the divergence in the expression for the Haniihonian density ()47p . we obtain 
the final expression for the Hamiltonian: 

i^ADM = j dxH^^,, = J dx{NH^ + N^H^). (49) 

We now rewrite the expression for Hx^ in terms of the irreducible components of 
momenta: 

H^ = \ P'-'P.u--, 7^, -y -eR. 

e \_ [n — l)[n — 2) 

In particular, the quadratic part of the momenta in Hj^ for n > 3 is consequently not 
positive definite. 

6 Secondary constraints 

To finish constructing the Hamiltonian formalism, we must analyze the consistency of pri- 
mary constraints ()43|) with the equations of motion. The phase space of general relativity 
is described by n{n+ 1) conjugate coordinates and momenta: (A^, p-*-), {N^,p^), {g^u^p'^^) 
on which the canonical equal-time Poisson bracket 

[N, p'^] = 6, [N^, j9'1 = 6;6, [9,.,p'n = 5P;j, (50) 

is given, where primed field variables are considered at a point x' = (x'^, . . . ,x'"~^). All 
fields are considered at the same instant t = x^. For brevity, we use the notation 



n— 1 „/n—l\ 



6 = (5(""^) {x - x') = 6{x^ -x'^)... 5{x- ' - X 

1 (51) 

for the {n — l)-dimensional 5- function and the symmetric combination of the Kronecker 
symbols in the right-hand sides of the Poisson brackets. We write all (5-functions on the 
right to distinguish them from a field variation. 

We now consider the Hamiltonian equations of motion for primary constraints PHj) : 

p^ = [p^, H^,^] = -H^, p'^ = [p^\ if^oM] = -H^. 

The consistency of the primary constraints with the equations of motion p-^ = 0, p^ = 
leads to the secondary constraints 

H^ = 0, H^ = 0, (52) 

where H± = H-^ and i^^ = g^^H'^. We note that the secondary constraints are not 
tensors but tensor densities of degree —1. Moreover, it is more convenient to consider the 
equivalent set of constraints with lowered index if^ instead of the constraints if ^. Below, 
we show that these constraints define the generators of coordinate transformations on the 
sections x^ = const and satisfy a simpler algebra. 

The constraints if^ are linear in the momenta and metric. The constraint H± is 
quadratic in the momenta and nonpolynomial in the metric g^i, because it depends on 
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the square root of the metric e and the inverse metric g^^'^ . The last circumstance raises 
essential technical difficulties in the perturbation theory. 

The secondary constraints are independent on the canonical variables {N,p^) and 
{Nfj_,p^), and they can be eliminated by considering the 77,(77, — l)-dimensional phase space 
of the variables Qf^iy and p'^'^ on which constraints (jH^ are imposed. In this case, the lapse 
and shift functions N and iV^ are regarded as Lagrange multipliers in the problem for the 
conditional extremum for the action 



S = jd-xip'^'^g^^-n. 



Because Hamiltonian ()49|) in general relativity is equal to a linear combination of the 
secondary constraints, we must compute the Poisson brackets of the constraints between 
themselves to analyze the consistency of secondary constraints (jH21) with the equations of 
motion. The constraint algebra in general relativity is well known, 

[H^, H'^] = -{H.g^'^ + H'^g'n^u, (53) 

[H^,H'^] = -H'J,, (54) 

[H^,H'J\ = -H,5,-H'^5,, (55) 

where we use the shorthand notation for the derivative of the (5-function, 

^. = ^^(^'-^)- 

Straightforward calculations of constraint algebra (j3^ - (j^ are very tedious. This 
algebra was first written by Dirac ^3] using symmetry considerations. The assumptions 
on the form of the constraints made in the course of derivation are not satisfied in general 
relativity, and the existence of the corresponding canonical transformation is now ques- 
tionable. Therefore, Dirac's derivation of the constraints algebra cannot be considered 
satisfactory. 

Two Poisson brackets (jH3) and (j33j) can in fact be found without straightforward 
calculations. For this, we consider the functional 

Tu = - fdxu'^H^, 

where u^{x) is an infinitesimal vector field. Calculating the Poisson brackets of the phase- 
space coordinates gfj_i, and p'^^ with T„ yields 

Kg^lu = [gfiu, Tu] = -d^ui^gpy - d^u'^g^p - u^dpgf,,,, 

This means that the functional T„, which is defined by the constraints if^, is the generator 
of general coordinate transformations on the hypersurfaces x° = const. We recall that 
the momenta p^^^ are tensor densities of degree —1. The algebra of general coordinate 
transformations is well known and defined by the Poisson bracket ()55p. We also can avoid 
computing the Poisson bracket ()54|) explicitly. Its form follows from the fact that the 
constraint H± is the scalar density of weight —1. Therefore, only the Poisson bracket 
(|^ must be computed. These calculations, being very cumbersome, were apparently 
first performed much later by DeWitt ^1. In the next section, we compute this Poisson 
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bracket after the canonical transformation, which casts the constraints into a polynomial 
form,essentially simplifying the calculations. 

We say that the constraints H^ are kinematical because they define only space diffeo- 
morphisms. They are also independent of the coupling constants in the action if there are 
any. The constraint H± is said to be dynamical because it governs the evolution of the 
initial data in time and depends essentially on the original action, in particular, on the 
coupling constants. 

For comparison, we write the Poisson brackets of the constraints H^ = with a 
contravariant index which are equivalent to the constraints H^^ = 0: 

We see that this seems more complicated than bracket (|33|). 

7 The canonical transformation 

The idea of the canonical transformation is as follows. The momenta p^^ are reducible and 
decompose into the traceless part and the trace PH|1 . Usually, working with irreducible 
components is more convenient for calculations because many terms automatically cancel. 
We pose the question: "Is it possible to perform a canonical transformation such that the 
irreducible components p^^ and p become new canonical momenta?" This question is 
nontrivial because the decomposition of the momenta involves the metric, whose com- 
ponents themselves are coordinates of the phase space. The answer to this question is 
negative because the Poisson brackets between the momenta are nonzero. For example, 
\P^^ iP'] 7^ 0- Nevertheless, there is a canonical transformation such that the new momenta 
are proportional to the irreducible components p^^ and p. Constructing this canonical 
transformation is our subject in this section. 
We consider the canonical transformation 

{9,u,pn - {k,.,pn,{p,p), (56) 

to the new pairs of canonically conjugate coordinates and momenta with the additional 
constraints on the coordinates fc^i, = k,yfj_ and conjugate momenta P'^'^ = P^^ 

|detV| = l, P'^^V = 0, (57) 

and p > 0. We choose the space integral as the generator of the canonical transformation 



I dx p'^k^^p^"' , meR, m^^O, (5^ 



depending on the new coordinates p-,k^y, old momenta p^^ , and the real parameter m. 
The old coordinates and new momenta are then given by the variational derivatives (see, 

e.g., na) 



'v-t^} (59) 



9i^v = 


^rp^lv P 


pfiu _ 




p = 


6F m 



(60) 
(61) 
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Computing the variational derivative with respect to k^^, we use the condition | det k^^\ = 

1, which restricts the variations k'^'^Sk^i, = 0, where k^'^ is the tensor density inverse to 

kfiu'- k^'^kya = 6^. Hence, the vanishing of momentum traces (j37jl follows automatically 

from the unit determinant condition for the density fc^i, for generating functional (j^Hj) . In 

Eq. (jnH), we use relation (jSHI)- 

In essence, the determinant of the metric raised to some power is singled out from the 

metric, 

1 

as a consequence of ()59p . For brevity in what follows, we also call the symmetrical tensor 
density with the unit determinant k^^, the metric. 

Variables ()5(i|l for n = 4 and m = 1/2 were considered in ^2], although the canonical 
transformation was not noted there. 

Straightforward calculations yield the expression for the scalar curvature of the section 
x° = const in the new coordinates: 

R = ^-—2 p2^(fc) + m{n - 2)pd^{k^''d,p) + m{n - 2) (^^^^ - l) k^"d^pd,p 

"(62) 
The "scalar curvature" for the metric k^^ takes a particular simple form, 

i?^'^) = d^k^^'' + h^'d.k.^d.k'"' - ^k'^'^d.k^^d^kP''. (63) 

This expression is not a scalar with respect to coordinate transformations of x'^, because 
kfj,i, is not a tensor but tensor density. But we note that the group of diffeomorphisms of 
sections x^ = const has a subgroup consisting of coordinates transformations of x^ with 
a unit determinant. The density k^,y is a tensor and R^''^ is a scalar with respect to this 
subgroup. 

As a consequence of the unit determinant of the metric, we obtain the components of 
the inverse metric k'^'^ as polynomials of degree n — 2 in the components k^iy, 

UfJ-^ — ^fJ.pi...p„-2 ±uai...<j„-2 I. I. 

where |5Mi--/^n-i| = x jg the totally antisymmetric tensor density of rank n — 1. Therefore, 
the scalar curvature R^''^ is polynomial in the metric fc^,y as well as in its inverse k'^'^. 
The dynamical constraint in the new variables becomes 



Tn.(n — 1) 

H± = p ^~ 



ppv p '_ p2 

^" m2(n-l)(n-2) 



p 2-^-'"-2 



p2i?(^0 + m{n - 2)pdf,{kf"'d^p) + m{n - 2) im- — ^ - 1 j k^'^'d^pd^p 



where P^^ = k^pk^^Pf". 

We now analyze the possibility of choosing the constant m such that the dynamical 
constraint becomes polynomial. Both expression in square brackets are polynomial in all 
dynamical variables. Because n > 3, we need the inequality m < to ensure a positive 
power of the density p before the first square bracket. In this case, the power of p before 
the second square bracket is negative. We therefore cannot ensure that the constraint H± 
itself is polynomial by choosing the constant m. But a constraint can be multiplied by 



an arbitrary nonzero factor without changing the surface defined by the constraint in the 
phase space. The power of p by which Hj^ be multiphed is minimum when the powers of 
p before the square brackets are equal We consequently obtain th equahty 

2 

m = . 

n-2 

Then multiplying the dynamical constraint 

K^ = p^H^ = eH±, (64) 

we obtain the equivalent polynomial constraint 

n-2 



K^ = pf^^P - ;^ ~ p'P' - p'R = 0, (65) 

4(n — 1) 

where the scalar curvature density 



p n — 2p2 n — 1 

is introduced. We note that the "scalar" curvature R^^"^ constructed for the metric density 
kfj,i, is not a scalar density. Therefore, using R instead of i?*-'^-' simplifies many formulas 
and calculations. 

Multiplying the dynamical constraint H± by the nonzero factor leads to modifying the 
Lagrange multiplier (the lapse function) by the inverse factor, 



n-l 



N^N = Np'^^^. (66) 

In turn, the modification of the Lagrange multipliers in general relativity is equivalent to 
coordinate changes that do not change the physical content of the theory. 

The kinematical constraints in the new dynamical variables remain polynomial: 

n-2, 
n — 1 



H, = -2V,{P\) - —^V,{Pp) = 0, (67) 



where V^i^ is the covariant derivative constructed for the metric 

2 

and indices are lowered using the tensor density P'^ ^ = P'^^kp^. It can be easily verified 
that 

V^p = 0, W^k.p = 0. 

Lowering and raising indices with the metric density k^u therefore commutes with the 
covariant differentiation operation. 

We note that the covariant derivative of a tensor density of degree deg = r in our 
notation is given by the expression: 

v^0 = 9^0 + rr^0, r^ = r,/ = !i^^. 

n — 2 p 
All new canonical variables are tensor densities of the degrees 

degkf,i, = -, degp = -, 

n — 1 n — 1 

n + 1 1 

degP^"' = ^, degP 



n — L n 
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and this should be taken into account for covariant differentiation. 

We now compute the basic Poisson brackets for the new canonical variables, which 
follow from explicit expressions (j3^ - ()61|) Only three brackets are nonzero: 

[p,P']=5, (68) 

Ku. p'n = {s^,: - ;7^ v^"'^) ^, (69) 

[P^^ P'p"] = —^{pt'^kP'' - pp''k^'')5. (70) 

Poisson brackets ()fJ9j) and (f7n|l do not have the canonical form for the phase variables. 
This occures because the fields k^^ and P^" are subjected to additional constraints (fKTjl . 
Because the Hamiltonian 



H= f dxiNK^ + N^H^) 



is polynomial in the new variables, the equations of motion are also polynomial. Straight- 
forward calculations with Poisson brackets (|F)H|) - (|7n|) yield the equations of motion: 

P = - T^P^.Np^P + ^pV.iV^ 
2[n — 1) n — 1 

ri — 2 ~ ~ ~ 1 

P =-, r-iVpp2 + 2NpR + 2pUN + -PV^N^' + N^V ^P, 

2{n — 1) n — 1 

k^, =2NP^, + V^N, + V,N^ - ——^k^^VpNP, 

p/^^ = - 2NP^'PP\ - p2 [Nk^Pk^^R,, —NRkP" + V^V^iV —k'^'aN] + 

\ n — \ n — 1 J 

+ VpiNPP^''') + ^—P^'^^VpA^^ - PPPVpN" - P^PVpN^, 
n — 1 



where 
and 



N^ = N'^k^,, W = kP-V,, □ = fc^'^V^V,, 



R = rW ^ ^ ~ 3 ^t^'^P _ jn - l){n - 3) d^pd^p _ n-3 r^^^d^p 
'"' '"'^n-2 p (n-2)2 p^ n-2 p" p 

{n-2f p2 '^n-2^'' p 

We now consider the constraints algebra. It is changed because we introduced the new 
constraint K±^ instead of the dynamical constraint H±. Simple calculations yield 

[K^, K'^] = -ip'H.k^'' + p"H'^k'nSu, (71) 

[K^,H'^] = -{K^ + K'^)6„ (72) 

[H,,Hl] = -HJ,-H'^6,. (73) 

Changes occur in Poisson brackets ()71|) and (f72|) as compared with the original algebra 
(J53|) - (j54j) . The second Poisson bracket, bracket (j72|) . has a kinematical origin and is 
defined by the fact that the new constraint is not a scalar function but a tensor density 
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of degree degK± = —2. Poisson bracket (fTTj) from direct calculations. We note that 
calculating that bracket is much simpler in the new variables than in the original ones. 

So far we have considered the metric density k^jy and its conjugate momenta P'^'^ 
together with additional constraints (j37|) . This is possible in the classic theory, but 
problems arise in the quantum theory of gravity. In the functional integral, considered 
in the next section, we integrate over all values of /c^jy and P'^'^. In principle, we can 
solve the constraints explicitly, but doing so is not interesting, because we then no longer 
have polynomials. Therefore, we describe the manifold N given by the coordinates kfj^^ 
and P'^'^ in detail. For simplicity, we assume that the coordinates take all possible real 
values, and the manifold N is consequently topologically trivial and diffeomorphic to the 
Euclidean space ]R"("~^\ We have the coordinate Poisson brackets (|69|) and (|7Up on that 
manifold, thus defining the Poisson structure. It can be easily verified that this structure 
is degenerate. This means that the manifold N is not a symplectic and is only a Poisson 
manifold (see, i.e., pTO]). Because the rank of the Poisson structure is n{n — 1) — 2, there 
are two functionally independent Casimir functions on the Poisson manifold N, 

Ci = det k^,, C2 = ^^P^^K^ (74) 

(introducing the constant factor l/(n— 1) in C2 simplifies several formulas in what follows.) 
Indeed, the Poisson brackets of the functions Ci and C2 with all coordinates are zero, 

[C,,2,k'^,] = [Ci,2,P'n = [Cl,2,p'] = [Cl,2,P'] = 

as a consequence of the definition of Poisson structure (jEHI), f^ - The Poisson brackets of 
these functions with an arbitrary differentiable function / e C^(N) then vanish, [Ci^2, f] = 
0, and Ci and C2 are therefore Casimir functions. The Poisson structure projected on 
sections V C N defined by the equations Ci^2 = const is nondegenerate. These sections 
are hence symplectic. 

It is always possible to choose local coordinates on the Poisson manifold N that are 
connected with the symplectic leaves Ci,2 = const. We let {qA,P^), A, B, . . . = 1, . . . , n{n — 
l)/2 — 1, denote the coordinates on these leaves. We choose coordinates ^a and p* such 
that the representation 

k,. = |Ci|^^^„ P'^^ = h^ + C2P^ (75) 

o o 

is satisfied, where the respective matrix elements k^jy{q) and P^^{p) depend only on g^ 
and p*. Obviously, such representations always exist. We choose the Casimir functions 
themselves as the lacking coordinates on N. We thus obtain the local coordinate system 

(fc^^Pn -- (gA,Cl),(p^C2). (76) 

Constraints (|F7J) have a simple form in the new coordinates, 

Ci = ±1, C2 = 0. (77) 

We choose the plus or minus sign for Ci if the space has the respective even or odd 
dimensionality. 

We consider the canonical transformation (j3^ - (jUT|) from a different standpoint. Strictly 
speaking, a canonical transformation considered in this section is canonical only between 
coordinates 

{9,..pn ^ {p,P),iq,,p^). (78) 
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The new phase space of general relativity in the considered case is the manifold M+ x M x V, 
where p G M+, P eM., and the submanifold V C N is defined by two values of Casimir 
functions (f77j) . The Poisson brackets on V have the canonical form by construction, 

The constraints are made polynomial by extending the space V to the Poisson manifold N 
with Poisson brackets (jUHjl - lfTOj) . When additional constraints (jHZI) are solved explicitly, 
the polynomiality is lost. This is not surprising. For example, electrodynamics contains 
constraints whose explicit solution even leads to a nonlocal action for physical degrees of 
freedom (see, i.e., jT7]). 

The Poisson brackets of the Casimir functions between themselves equal zero [Ci, C2] = 
0. From the standpoint of the Hamiltonian formalism, they could be regarded as the first- 
class constraints generating gauge transformations. But these transformations are trivial 
because the Poisson brackets of Casimir functions with all phase-space coordinates vanish. 
This is possible only on a Poisson manifold with a degenerate Poisson structure. There 
are no Casimir functions on a symplectic manifold. 

The Poisson manifold N can be equipped with the second, now canonical, Poisson 
bracket. With respect to this new canonical Poisson structure, the submanifold V is de- 
fined by two second-class constraints (j77|) . Then the original degenerate Poisson structure 
(I69|) . (j7(jp is just the Dirac bracket with respect to the canonical Poisson structure on N. 

We also note the following. In gravity, we assume that a space-time metric gap has the 
Lorentzian signature and therefore is nondegenerate. In quantum gravity, the functional 
integral is integrated over all independent metric components, and this property can be 
taken into account only by restricting the integration domain. When the phase space 
is extended to a Poisson manifold, the integration domain is extended to the Euclidean 
space. The nondegeneracy of the metric is automatically provided by the presence of 
5-f unctions in the integrand. 

8 Generating functional for the Green functions 

Here, we summarize the results calculating in the most interesting case of the four di- 
mensional space-time and write the explicit expression for the generating functional for 
Green's functions. The polynomial Hamiltonian formulation of general relativity on the 
Poisson manifold M+ x M x N with coordinates (p, P) e M+ x M and {k^^, P'"') e N is 
given in the preceeding section. The dimensionality of this manifold for n = 4 equals 14. 
The Poisson structure is defined by the nonzero Poisson brackets 

[p, n = s, 

[kfiu, P \ ~ \ ^tiy ~ q ^Mi^^ ) "' 
o 

It is degenerate and has rank 12, which coincides with the dimensionality of the phase 
space of general relativity. There are two Casimir functions (I74j) on N. The section 
Ci = — 1, C2 = is a symplectic submanifold V G N and defines the phase space of 
general relativity. 
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The action of general relativity in the new coordinates has the form 

^HE = I dx{Pp + P^'k^, -H- d^B^), (79) 

where the Hamiltonian density H is equal to a linear combinations of constraints, 

H = NK^ + N^H^, 

and N and N^ are Lagrange multipliers. To action (j79j) . we added the boundary term 
d^j^B^^ on a space section x° = const, which is written as the divergence of some function 
of the canonical variables B^{p,k^y,P,P^"). We briefly discuss the necessity of adding 
this important term in the action below without specifying its form. The constraints 

K^ = P^^P,, - ^p'P^ - p^R^'^ - 2pd,{k>^''d,p) + hf^'d^pd^p, (80) 

H^ = -2d,{P'"^k,^) + P'^'d^k,, - h,{Pp) + Pd,p, (81) 

are polynomial first-class constraints and satisfy the algebra (J71|) - (j73p . The scalar cur- 
vature w-^^ for the metric density k^u with unit determinant has the form ()63|). The 
constraint if _l is quadratic in the momenta and the variable p. It is a fifth-order polyno- 
mial in metric density k^y (and its partial derivatives). The constraint if^ is linear in the 
momenta and also in the coordinates. 

The expression for the generating functional for Green's functions as a functional 
integral over the phase space ^H] is easily generalized on a Poisson manifold. For brevity, 
we introduce a new notation for the secondary constraints and Lagrange multipliers: 

{Ha} = {K^, H,}, {iV-} = {N, m}, a = 0, 1, 2, 3. 

We now fix the invariance under general coordinate transformations using four gauge 
conditions F"" = 0. The gauge is assumed to be canonical, 

det[Ha,F']^0. 

The canonical form of the generating functional for the Green's functions for the metric 
is given by the functional integral up to a normalization factor 13j: 

Z( J) = J Dig,,)DipnD{N'') ^^P U j dx{p^''g,, - N'^H, - d^B^^ + g.^J^ } x 

X det[ii„,F^]J]5(iJ,)n^(^"), (82) 

a a 

where 

X, fj,<u X, fi<u X, a 

and the J'^'^ are the sources for the metric. The functional integral Z{J) is the generating 
functional only for "coordinate" Green's functions because the sources are written only 
for the space metric |T7j. Because the Jacobian of any canonical transformation and. 
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in particular, of transformation (j78p is equal to unity, the expression for the functional 
integral can be rewritten in the equivalent form 

Z{J) = J Dip)DiP)Diq,)Dip-)D{N'^)x 

X exp I ^ /" dx{Pp + p^q^ - N^Ha - 5^5^ + pJ + q^J"" ) > x 



X 



6et[Ha,F']\{5{H,)\{5{F'^) 



The constraints Ha are nonpolynomial in coordinates q^ and momenta p^ in this 
form. The integration must be extended over the whole Poisson manifold N to make the 
constraints polynomial. We provide this by introducing two additional (5-functions: 

Z{J)= j D{p)D{P)D{k^,)D{P^'^)D{N'')x 

X exp I ^ /" dx{Pp + P^'^'k^, - N'^Ha - df,B^ + pJ + k^J^"") | x (83) 



det [Ha,F']6{Ci + l)6{C2)l[6{Ha)l[6{F^). 



Two new (5-functions remove the integration over the additional variables and restrict the 
integration over the Poisson manifold N to integration over symplectic section ()77|) . To 
perform the corresponding integration, we must perform coordinates transformation (j76|) . 
We must prove that the Jacobian of the coordinates transformation (j76|) 



D{q,)D{C,)D{p')D{C2) = D{k^,)D{P 



flU\ 



d{q,,Cup\C2) 



d{k^u, P'^") 



equals unity on constraint surface (f77|) to prove that the proposed expression for generating 
functional (jHHj) over the Poisson manifold is equivalent to the original functional integral 
(j82j) over the phase space. This can be easily done. Let the coordinates ^a parametrize the 
matrix k^^, with a unit determinant arbitrarily. Then 



dik^.,P 



'fJ.U] 



= 




d{p\C2) 
d{Pt^^) 



because the matrix 



9(P^^) 



0. 



As a consequence of the definition of canonical momenta (jUOj) . we obtain 

. 5F ^ 5F dk^, 

6qA 6k^u dqA 



P 



rik 



dqA 



Representation (f7^ for the metric fc^i, yields 



"/lU 



Therefore, 



dC2 

Qpfiu 



{n - l)Ci 



u 

n — 1 



dk 



/lU 



dCi 
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Consequently, 



d{p\C2) dk,. 



This implies that for the constraint surface, the modulus of the Jacobian of the coordinates 
transformation is equal to unity: 



d{q,,C,,p\C2) 



9(A;^.,P^-) 



Ci=-1 



This motivates introducing the numerical factor in the Casimir function C2 in Eq. ()74|) . 
Without it, the Jacobian of the coordinate transformation would equal some nonzero 
constant, which could be included in the definition of the normalization factor of the 
generating functional. 

We now say a few words about including the boundary term d^B^ in the action of 
general relativity. The original expression for the generating functional (jH^ is justified 
by assuming that the functional integral over just physical degrees of freedom with the 
unit measure appears after solution of all the constraints and gauge conditions. The 
Hamiltonian on the constraint surface becomes zero in this case, N"'Ha = 0. On the 
other hand, we know that the dynamics of the physical degrees of freedom are nontrivial. 
A possible way out from this contradiction was proposed in 0. If the gauge condition 
depends explicitly on time, then the nontrivial Hamiltonian for the physical degrees of 
freedom on the constraint surface arises from the kinetic term p'^'^g^u- But this is not a 
unique possibility. To obtain a nontrivial Hamiltonian for the physical degrees of freedom 
in the canonical gauge, which does not depend on time explicitly, the boundary term was 
added to the action J7j. The importance of the boundary term in general relativity is 
now universally recognized, but its role still remains obscure in many cases because of 
large technical difficulties. In two-dimensional gravity where the constraints are solved 
explicitly, it was proved that the Hamiltonian for the physical degrees of freedom appears 
from the boundary term for gauge conditions, which do not depend explicitly PE|- Namely, 
the term d^B^^ on the constraint surface is equal to the Hamiltonian density for the 
physical degrees of freedom and is not equal to the divergence of some function. This 
statement is local and is independent of whether the universe is closed. 

We note one more important point. For an asymptotically flat space-time, the volume 
integral of <9^-B^ is equal to the surface integral of B'^ and coincides with the mass of the 
Schwarzschild solution. That is why this integral was proposed as the definition of the 
total energy of a gravitational field |2]. It is usually assumed that the total energy of a 
closed universe vanishes because it has no boundary, and no surface integral arises. This 
is true if the constraints admit smooth solutions on compact manifolds for the nonphysical 
degrees of freedom in the general case. The example of two-dimensional gravity (which 
includes spherically symmetric solutions of general relativity) shows that constraints do 
not admit smooth solutions on a circle in general. In this case, we must make a cut on a 
compact space and add a boundary term there in order to pose the variational problem. 
This yields a nontrivial expression for the total energy of a closed universes. We can say 
this in other words. The divergence d^B^ on the constraint surface expressed in terms 
of the physical degrees of freedom is no longer the divergence of some function. This 
Hamiltonian density for the physical degrees of freedom has the same form independently 
of whether the universe is closed or not. Therefore, this expression can be taken as the 
definition of the energy density of a gravitational field. Of course, this definition is not 
covariant and depends on the choice of a coordinate system. 
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9 Homogeneous and isotropic universe 

We consider the Friedman universe [fl^ 120] to demonstrate the new variables introduced 
in Sec. |7| First, we recall the derivation of the equations in the Lagrangian formalism. 
For comparison, we then reformulate the model in the Hamiltonian language in the old 
and new variables in which the Hamiltonian becomes polynomial. 

In the Friedman model of the universe, which is the basis for most contemporary 
cosmological models, we assume that the space is a Riemannian manifold of constant 
curvature at each instant (all sections x° = t = const). This assumption corresponds to 
a homogeneous isotropic universe. The metric satisfying this requirement is 

ds'^ = df + a^g^^dz'^dx", (85) 

where a = a(t) is the scale factor depending only on time, and g^^, is a (negative-definite) 
constant-curvature space metric, which is independent of time by assumption. The specific 
form of the constant-curvature space metric g^^ depends on the coordinate system and is 
not important for the following consideration. For simplicity, we restrict ourselves to the 
most important case of four-dimensional space-time (n = 4). 
The metric of the space-time satisfies Einstein's equations 

where we introduce the cosmological constant A and the matter energy-momentum tensor 
Tq,^. The energy- momentum tensor in the comoving frame has the form J2I] 
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where £ and V are the respective energy density and pressure of matter. For a homoge- 
neous isotropic universe in the chosen coordinate system, these densities depend only on 
time 8 = S{t) and V = V{t). 

We also assume that the matter equation of state is given by some dependence of 
energy density on pressure, 

£ = £{V). (88) 

Because energy-momentum tensor (|87|) was not obtained by varying some invariant ac- 
tion for matter fields with respect to a metric, the energy density £ cannot be an arbitrary 
function. Indeed, the covariant divergence of the left hand side of Einstein's equations is 
identically zero as a consequence of the Bianchi identities. Therefore, Einstein's equations 
yield the equation for the energy-momentum tensor: 

V,T-^ = 0. 

For metric ()85p and energy- momentum tensor ()87|) . these four relations reduce to one 
nontrivial equation 

£ + —{£ + V) = Q. (89) 

a 

For a given equation of state (jHHj) . we have one first-order differential equation, which we 



rewrite in the form 



d£ da , ^ 

3—. (90) 



£ + V{£) a 
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Solving this equation yields the energy density S as a function of the scale factor a. 

We obtain the equation for the scale factor because metric (jH^ must satisfy Einstein's 
equations ()86p . Simple calculations yield the expressions for the Einstein tensor 

1„ ^Ko-d^ 



Rq R = 3 



2"~ " a2 ' 
Ro^ = 0, 

R/ - hibl = -^ {2ad + d'- K,) 5;, 

where Kq in our notation is the scalar curvature of a three-dimensional sphere {Kq = 1), 
Euclidean space {Kq = 0), or one-sheet hyperboloid {Kq = —1). We can now easily see 
that Einstein's equations lead to only two nontrivial equations on the scale factor: 

3^^^ + A + -^ = 0, (91) 

-^(2ad + d^ ~Ko)+A--V = 0. (92) 

a^ 2 

It is easy to verify that Eq. (jU^ is a consequence of Eqs. (PT|) and (|H^ because Einstein's 
equations are linearly dependent as soon as Eq. (jH^ is satisfied. Therefore, Eq. (j^ can 
be dropped, but we do not do this, because it is needed for the canonical treatment of 
the Friedman universe. 

We thus find the dependence £ = £{a) of the energy density on the scale factor by 
solving Eq. (jnUjl for a given equation of state (fHH|l . Substituting this function in Eq. 
()91|). we obtain the first-order ordinary differential equation for the scale factor. This 
is precisely the main equation in the standard cosmological models for a homogeneous 
isotropic universe. 

We start with the canonical formulation. For metric (|85|) . we have 

iV = 1, N^ = Q, g^^ = a^g^^, 
External curvature tensor (|32|) and the volume element are 

K = -adg K = -3-. 

^ a 



a e, e = y | det^f 



y fiui 



Canonical momenta ()44j) conjugate to the space metric become 

p^" = -2deg . 

The traceless part of these momenta is identically zero, 

P^'•' = 0, p= -Ga^de. (93) 

The Hamiltonian density is given by the single dynamical constraint ()48j) because the 
shift function for the Friedman universe is equal to zero 



Hi =-6e 



(d^ - Ko)a - -Aa^ - -£a^ 
3 6 
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(94) 



where we take the contribution from cosmological constant and matter fields into account. 
We cannot insert the momentum trace p instead of the time derivative a in this expression, 
because it is not the variable conjugate to the scale factor. 

To find the momentum conjugate to a, we rewrite Lagrangian (jl^ for metric ()85p . 



•i-'AriM — 06 



6 D 



This expression can be integrated over the space because all the dependence on the space 
coordinates is contained in the volume element e. Dropping the constant factor —6V, 
where V is the volume of space (infinite for the Euclidean space and one sheet hyper- 
boloid), we obtain the Lagrangian for the scale factor 

L = {a^ + Ko)a + -Aa^ + -Sa^, (95) 

3 6 

which is independent on the space coordinates. This is the standard Lagrangian for a 
point particle moving in a one-dimensional space with the coordinate a G M^. 

Deriving the Lagrangian for the scale factor, we dropped the negative factor to obtain 
the positive sign of the kinetic term aa^. We thus changed the total sign of the action 
and hence the sign of the energy. We note that the contribution of the kinetic term of 
the scale factor to the energy is negative. 

The expression for the momentum conjugate to the scale factor follows from La- 
grangian (jHS)): 

dL ^ . 

Pa = ^^ = ^aa, 
oa 

which differs from momentum trace (j^Hj) by a factor and is independent on the space 
coordinates. The Hamiltonian for the scale factor corresponding to Lagrangian (j95|) is 

H = ^pI- Koa - \Aa^ - |a^ (96) 

We see that this expression for the Hamiltonian coincides up to the factor —6V with 
the expression obtained by integrating dynamical constraint (|94p over the space. This 
observation is nontrivial because some equations of motion may be lost when specific 
expressions for field variables are inserted into the Lagrangian. 

The Hamiltonian equations for the scale factor are 

1 

2« ^ (97) 

Pa = -^pI + Ko + Aa^ + -Va\ 

where we use Eq. (jHIH) to calculate the Poisson bracket 

[Pa,t\ = [Pa,a\—- = -3 . 

da a 

It can be easily verified that Hamiltonian equations (jUTj) are equivalent to the second- 
order Lagrangian equation (j^^ . and Hamiltonian (j^Uj) is proportional to the left-hand 
side of Eq. dHU). 
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We have thus formulated equations for the Friedman universe in the Hamiltonian lan- 
guage. In contrast to the Hamiltonian particle dynamics, we have an additional constraint 
along with canonical equations of motion ()97|) . 

H{a,pa) = 0. (98) 

In other words, we only seek those solutions of the equations of motion for which the energy 
is equal to zero. This problem is self-consistent because the energy is conserved. We have 
thus formulated equations ()91|1 and (jH^ for the scale factor in the Hamiltonian form and 
proved that one constraint (j98|) on the scale factor and the corresponding momentum is 
lost if the expression for the metric (|H3j) is substituted not in Einstein's equations but in 
the action. 

Hamiltonian (jUU]) and equations of motion ()97|) are nonpolynomial in the scale factor. 
We show what happens with the equations under the canonical transformation described 
in sectional As a consequence of Eqs. (jKn|) - (jHT|l . we obtain expressions for the canonical 
variables 

(99) 
(100) 

external curvature 



and volume element 



after the canonical transformation. 

To simplify calculations, it is easiest to separate variables by extracting the factor q{t) 
depending only on time from p, 

o2/3 

p{t,x'')=q{t)e . 

The Lagrangian density can then be integrated over the space as previously, and the 
Hamiltonian reformulation of the equations reduces to the redefinition of the scale factor 

q = a^. (101) 

Lagrangian (j^3j) for the new variable is 

4:y/q 3 6 

The momentum conjugate to the new dynamical variable q{t) is 

The corresponding Hamiltonian contains a nonpolynomial factor. 



p = 


a'p/', 


f^fiu Vl-ii^^ 


p = 


Vp 


pM^ = 0, 


K^^ = 


1 


2p 




e = p 


3/2 



H = ^{pl-Ko--Aq- -Sq 
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During the construction of the polynomial Hamiltonian formulation, the dynamical con- 
straint was multiplied by factor ()64|) . As a result, we obtain a new constraint, which is 
now polynomial, and the new Hamiltonian 



K = q"'H = q' {pI -Ko- Uq - Isq^ 



(102) 



To preserve the Hamiltonian form of the equations of motion we must also redefine time 
as well t —>■ T, where the new parameter is defined by the differential equation 

dt ^ 

Redefining time corresponds to redefining Lagrange multiplier (jM]) . The equations of 
motion for Hamiltonian ()102|) are 

d(l o 2 

dr ^ \3 12 4 y g ' 

where we once again used Eq. (jHOI)- The second term in the second equation can be 
dropped because the model contains the constraint 

K{q,pg)=0. (103) 

Thus in the new variables, we have a "particle" described by the coordinate g(r) G M_(_ 
and momentum Pg(r) G M with Hamiltonian ()102j) and with constraint ()103p imposed on 
the canonical variables. The Hamiltonian and equations of motion are polynomial (if the 
energy density S{q) is polynomial in q) and equivalent to the original Einstein's equations 
for the scale factor (jHH), (El- 

For the homogeneous isotropic universe we can go even further and eliminate the 
common factor g^ from Hamiltonian ()102|) by redefining time t —>■ t', where 

dr' = q-^''^dt. 

We then have a "particle" with the simple Hamiltonian 

K' = pI-Ko- ^Ag - Isq. 



3 ' 6 



and the corresponding equations of motion 
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dr' 
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The space curvature Kq now does not contribute to the equations of motion at all. The 
system of equations of motion is solved under the zero constraint on the Hamiltonian, 
K' = 0. Just as for the original system of equations ()91|) and (j^J^ . the Hamiltonian 
equations of motion in the Lagrangian form follow from the equation K' = 0, where 
instead of the momentum, we must substitute its expression in terms of the time derivative 
of the scale factor (jl04|) . This equation is equivalent to Eq. (jHH). 
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For dust matter, we have V = 0, and Eq. (^01) can be easily integrated, 



c = -^, M = const. 



which provides a massless particle moving in the potential 



1 . M 

-An 

3 ^ 6g2 



Introducing new canonical variables thus allows reformulating the equations for the 
Friedman universe in polynomial form. Because Eq. ()9H1 has been well studied in recent 
decades, the new formulation is unlikely to yield new results for the analysis of classical 
solutions in the considered case. It may be useful for constructing a quantum model of the 
Friedman universe, but this question requires a separate investigation and is beyond the 
scope of this paper. Our consideration here is intended only to illustrate how the general 
method works in a simple case where all steps can be verified by simple calculations. 

10 Conclusion 

We have demonstrated that considering the determinant of the metric and the conjugate 
momentum as independent additional variables leads to a Hamiltonian formulation of 
general relativity in which all the constraints are polynomial. The model is formulated 
in polynomial form not in the phase space but on a Poisson manifold where the Poisson 
bracket is degenerate. We stress that the resulting model is equivalent to general relativity. 

In the new variables, the canonical momenta are proportional to the irreducible com- 
ponents of the momenta in the standard metric formulation. This property essentially 
simplifies the calculations, in particular, the calculations of the Poisson bracket of the 
dynamical constraint. 

The proposed canonical formulation of general relativity allows writing the functional 
integral on the Poisson manifold. We proved that this integral is equivalent to the func- 
tional integral over the phase space. The advantage of the new expression for the generat- 
ing functional for the Green's functions is that the action and all arguments of ^-functions 
are polynomial in independent variables. This leads to the presence of only a finite number 
of vertices in the diagram techniques. This seems to simplify calculations in the quantum 
theory of gravity. 

As an example of using new variables, we considered the Friedman model of the 
universe. In these variables, the Hamiltonian and the equations of motion for the scale 
factor are simplified and take a polynomial form. 

A similar transformation of variables in the configuration space leading to a polynomial 
Hilbert-Einstein action was proposed in p5? and recently rediscovered ^5] . 
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